Multipartite entanglement in quantum spin chains 
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We study the occurrence of multipartite entanglement in spin chains. We show that certain 
genuine multipartite entangled states, namely W states, can be obtained as ground states of simple 
XX type ferromagnetic spin chains in a transverse magnetic field, for any number of sites. Moreover, 
multipartite entanglement is proven to exist even at finite temperatures. A transition from a product 
state to a multipartite entangled state occurs when decreasing the magnetic field to a critical value. 
Adiabatic passage through this point can thus lead to the generation of multipartite entanglement. 
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Quantum entanglement is a valuable resource with po- 
tential applications including quantum frequency stan- 
dards 0|, quantum cryptography Q and quantum tele- 
portation |3( . Recent progress in experimental techniques 
allows to generate and to control multipartite entangle- 
ment 0, yj ■ This motivates studies of quantum entan- 
glement in more complex physical systems, for example 
chains of interacting spins. Here we show how passage 
through avoided crossings in spin chains can be used to 
generate multipartite entanglement. 

Quantum spin chains have been extensively studied in 
the context of quantum information science, in particu- 
lar their use as quantum wires and as simple quantum 
processors |?J . Bipartite entanglement in spin chains has 
been also thoroughly investigated H H E IH IH 13 
Il4lll5|. In contrast, very little is known about multipar- 
tite entanglement in such systems. A notable exception 
is work by Wang, who performed a numerical analysis of 
multipartite entanglement in the Heisenberg model |l6| . 
and by Stelmachovic and Buzek, who showed that in the 
limit of infinite coupling strength the ground state of the 
Ising chain is locally unitarily equivalent to the iV-partite 
generalization of the GHZ state [ijj . More recently San- 
tos showed that the presence of defects in a spin chain 
governed by the XXZ Hamiltonian could be exploited 
to generate bipartite and tripartite entangled states be- 
tween selected spins |l8|. 

Here we analyze a quantum spin chain with an arbi- 
trary number of sites and demonstrate the existence of 
multipartite entanglement of the W-type (see Eq. i|10|) 
below) in the ground state of the XX Hamiltonian with 
suitable transverse magnetic field. We find the critical 
value of the magnetic field at which the ground state 
undergoes a transition from a product state to this mul- 
tipartite entangled state. Our work therefore provides a 
method for generating a multipartite W-type entangled 
state by tuning a single external parameter, namely the 
global magnetic field, while the inter-spin interaction re- 



mains constant. We refer to spins as qubits and use |0) 
and |1) to denote spin-up and spin-down states, respec- 
tively. 

Let us start with a simple example, where the tran- 
sition from a product state to a multipartite entangled 
state can be easily seen. We consider the XXZ model 
for a chain of 3 spins in an external magnetic field and 
with periodic boundary conditions. The corresponding 
Hamiltonian is given by 
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(1) 



where J is the coupling in the x and y directions and 
J + A is the coupling in the z direction. The parame- 
ter A quantifies the anisotropy in the interaction; for the 
Heisenberg interaction A = 0. N.B. this model exhibits 
the same symmetry as the XX model with N sites which 
will be discussed later. In both cases the z-component of 
the total spin, of ot = J] i of , commutes with the Hamil- 
tonian H . Thus, the eigenstates of H are superpositions 
of states with a fixed number of up-spins. 

In Eq. Q the term proportional to A commutes 
with the terms proportional to J. Thus, the eigen- 
states of the Hamiltonian Hxxz coincide with those of 
the isotropic Heisenberg model. There are four non- 
degenerate eigenstates, given by llll),^) = (| 110) + 
|101) + |011»/V3, |W0 = (|001) + |010) + |100»/V3, 
and 1 000). The corresponding eigenvalues are Em = 
3 J + 3B + 3A, E w = 3J + B - A, E w = 3J - B - A, 
and -Eooo = 3J — ZB + 3A. The remaining eigen- 
states are doubly degenerate and are given by \WW) = 
(|110) +e 27 " fc / 3 |101) +e- 27rlfc / 3 |011))/ % /3 with k = 1,-1 
and eigenvalue £^(*o = J + B — A, and |W^) = 
(|001)+e 27rifc / 3 |010)+e- 2 " fc / 3 |100))/V3 with eigenvalue 
E W (k) = J — B — A. 

The relative values of the parameters J, B and A deter- 
mine which of these states is the ground state. Let us first 
consider the isotropic case, A = 0. When the value J is 
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FIG. 1: Role of the anisotropy: Energy levels for the ferro- 
magnetic XXZ model with anisotropy A as a function of the 
magnetic field B, for the case N = 3. Dashed lines are for 
A = and solid lines for A > 0. For B > B c the levels are 
ordered from bottom to top as Eooo,Ew, Eyy, £111 (see main 
text for the definition), both for A = and A > 0. 



negative, i.e. in the ferromagnetic case, the ground state 
is always a product state, regardless of the value of B. 
Changing the sign of the magnetic field B simply leads 
to the relabeling <-> 1 of each qubit. In the following, 
we will always consider B to be positive. For positive 
J, i.e. the antiferromagnetic case, and for B > J the 
ground state is also a product state. However, for B < J 
the ground state lies in the subspace spanned by iWfe), 
with k = 1,-1. This subspace contains states with W- 
type tripartite entanglement and biseparable states, i.e. 
states in which one qubit is unentangled. 

Let us now turn to the anisotropic case, A ^ 0. For 
negative values of J and negative A, the ground state is 
again given by a product state, regardless of the value 
of B. However, the case of negative J and positive A 
leads to interesting results. For high values of B, the 
spins are aligned, and the ground state is |000), i.e. a 
product state. When one decreases the value of B, at 
B c = 2 A the ground state changes to \W), namely a 
genuinely tripartite entangled state. The essential role 
of the non-vanishing anisotropy A for the existence of 
the ground state transition from a product state to a 
tripartite entangled one is illustrated in Fig. ^ 

The above considerations hold if the spin chain is at 
zero temperature. However, at finite but low tempera- 
tures, the state of this three-spin chain can be shown to 
be still genuinely tripartite entangled. To see this, one 
expands the density matrix of the system at a tempera- 
ture T, which is given by its Gibbs state, as follows 

p cx e-^ Eooo \000)(000\+e- pEw \W)(W\ 

+ e -' 3S ^|V?)(H>|+e- /3£lll |lll)<Hl| + ... , (2) 

where [3 — 1/kT and k is the Boltzmann constant. In the 



vicinity of B c = 2 A, the two lowest eigenvalues are -Eooo 
and E\y . Hence by retaining only the two leading terms 
of the expansion (J2J , which is a good approximation for 
sufficiently low temperatures, the state of the system can 
be approximated as 

Q = p\W){W\ + (l-p)\QQG){mQ\ , (3) 

where p = 1/(1 + exp[=F2/3 (2A - B)}) for B c < 2A and 
B c > 2A, respectively. This state was first studied in 
the context of approximate quantum cloning [l9| . It is a 
genuinely tripartite entangled mixture, which cannot be 
written as a mixture of biseparable states for any value of 
p =/= 0. To prove this, let us assume the contrary, namely 
that g is biseparable. A biseparable state g of a tripartite 
system ABC can be decomposed as [2(J 

8 = E (tfW + Pf^vf + Pir^ni) , (4) 

where pf , pf , p£ are probabilities with J2ijk(P? + pf + 
p£) = 1. In the first term on the RHS of Q, t^ b de- 
notes the joint density operator of the subsystems A and 
B, while r)f denotes the density operator of the subsys- 
tem C (and analogously for the two subsequent terms). 
The decomposition implies that there exists a bisepa- 
rable state of the form \tpf B ) \4>f) in the range of g. This 
in turn implies that there should exist non-zero coeffi- 
cients a and such that \ipf B )\<j>f ) = a\W) + /3|000). 
However, it is straightforward to see that there are no 
biseparable vectors in the subspace spanned by \W) and 
|000), apart from the vector 1 000) , which however corre- 
sponds to the trivial case a = 0. Hence one arrives at 
a contradiction. Let us already note here that this ar- 
gument can be generalised in a straightforward way to 

parties: any mixture of an N-party W state and the 
product state |00...0) is genuinely multipartite entangled. 

Let us now turn to the general case of a spin chain 
with N sites. Here we will focus on the choice A = — J, 
namely the XX model, with periodic boundary condi- 
tions, described by the Hamiltonian 

N 

H XX = J(af ® of +1 + at ® af +1 ) + B af) . (5) 
»=i 

Again, we consider the ferromagnetic case J < 0, and 
B > 0. As in the case of the XXZ model with three 
sites, [Hxx,a z to t] = holds, where of ot = £ili°f- 
Hence, the Hilbert space decomposes into invariant sub- 
spaces, each corresponding to a distinct eigenvalue of 
of ot which we denote as m. It represents the total num- 
ber of excitations (down-spins) in the chain. There are 
N + 1 such subspaces, corresponding to spin configura- 
tions with ?7i = 0,1,2,...,A. We refer to the subspace 
corresponding to a particular m as the m-excitation sub- 
space. The subspaces for which m — and m = N 
are one-dimensional. The eigenstates of the Hamiltonian 
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H in these subspaces are given by the product states 
|000...0) and |lll...l), respectively. The correspond- 
ing eigenvalues are 



E 



(o) _ 



-NB 



and E {N) = +NB 



(6) 



The energy eigenvalues for the m-excitation subspace can 
be obtained using the standard fermionization technique, 
introduced by E. Lieb et al and S. Katsura 21] . The 
eigenvalues are labeled by m distinct quantum numbers 
fci, k2, ■ ■ ■ , km € {1, 2, , . . . , TV} and can be expressed as 



E 



(m) 

fei > ■ "j^n 



4.7 



cos 



, 27rfci ■ 

"tv - 

, 27rfc„ 



cos( ^v~ ) + 



+ cos(— -) -(TV-2m)£. (7) 



TV 

In particular, for the single excitation subspace, 



(i) _ 



4Jcos(^)-(TV-2)T3, 



(8) 



where k = 1,...,TV. The corresponding energy eigen- 
states are given by 



1 N 

1 \ ^ 2-rrikn 
1^ 



/TV, 



(9) 



where the states |n), with n = 1, 2, . . . , TV, correspond to 
spin configurations in which all spins are up, apart from 
the spin at the site n which is down. For example, for 
TV = 3 we have |1) = |100), |2) = |010) and |3) = |001). 
These states form a complete basis in the single excitation 
subspace. 

Let us now look at the ground state of the chain as a 
function of the magnetic field B. If the chain is held in a 
very strong magnetic field B ^> 3 then its ground state 
is the product state of the form |000...0) and energy 
= —NB. However, when we decrease B, keeping 3 
fixed, then the lowest energy single excitation eigenstate 
\<Pn), which from now on we will label as |Wjv), 



\Wi 



1 



N / 



1 



|2) + |3) 



|TV> 



(10) 



becomes the new ground state with energy Ejp — 4 3 — 
(N—2)B. The crossover from the product state |000 ... 0) 
to the entangled state \ Wn) occurs at the critical value of 
the field B c = —2 J, which is independent of TV. The new 
ground state \Wn) is a generalization of the tripartite 
entangled VT"-state to TV spins and is genuinely TV-partite 
entangled. 

This is the first crossover because in our particular case 
of J < the lowest energy in the m-excitation subspace 
(for m > 2) is lower bounded by 4 3m — (TV — 2m)B (see 

Eq. @), which can be written as E$ +2(m- 1)(2 J+B) 
and is certainly greater than E$ for all values of B > 
-23 . 



Note that the value B c = —2 J corresponds to the criti- 
cal point for the Bose Hubbard model with infinite on-site 
repulsion energy of the bosons |22| . namely to the quan- 
tum phase transition from the Mott insulator to super- 
fluid phase. We also point out that, as explained below 
Eq. at finite temperatures and near B c = — 2J there 
remains genuine multipartite entanglement. 

In order to generate the state \Wn) from the state 
1 00 . . . 0) by lowering the magnetic field B from an initial 
large value, it is necessary to turn the level crossing at 
B c = —23 into an avoided crossing. An avoided crossing 
can be realized by adding a small perturbation to the 
Hamiltonian Hxx, e.g., a term B' J2iLi a ii w i tn B' << 
B. In the presence of an avoided crossing, a transition 
from 1 00 . . . 0) to \Wu) can be achieved by lowering the 
magnetic field slowly enough, so that the ground state 
adiabatically follows [23j . The required rate of change of 
B depends inversely on the gap of the avoided crossing. 
Using degenerate perturbation theory, this gap is found 
in first order to be equal to 2B'yN. 

Before concluding, let us make a brief remark about 
the ground state for values of the magnetic field in the 
range B < — 23 . It turns out that the ground state 
changes successively from one excitation number to the 
next-higher one when decreasing the value of the mag- 
netic field B. To show this, let us calculate the crossings 
of the energy levels for different excitation numbers (for 
large TV). It is straightforward to see that the lowest en- 
ergy in the m-excitation subspace, for odd values of m, 
is given by 



E {m) -41 



[W2] „ . 

1+ E 2cos (^r)] -{N-2m)B (11) 



whereas, for even values of m, it is given by 



E 



(m) _ 



4J 



to/2-1 . 



(N-2m)B. 

From this it follows that the crossing T5( m + 1 ) = oc- 
curs at B { c m) = -2J[1 - (2m 2 7r 2 )/TV 2 ]. Thus, with the 
addition of a suitable small perturbation, and by adiabat- 
ically decreasing the magnetic field, transitions through a 
cascade of ground states with increasing excitation num- 
bers can be obtained. 

In conclusion, we have shown that a multipartite en- 
tangled state of the W-type occurs naturally as a ground 
state in the ferromagnetic XX spin chain with TV sites in 
an external magnetic field. W-states are a useful resource 
for quantum information processing tasks, e.g. quan- 
tum teleportation and dense coding [241 ]. Our analysis 
suggests a new method of generating an TV-partite en- 
tangled W-state, namely by driving the chain adiabat- 
ically through an avoided level crossing. This amounts 
to preparation of the initial product state |000 ... 0) in a 
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strong magnetic field, B 3> J and then slowly reducing 
the strength of the field, in the presence of a small per- 
turbation, until it passes through the value B c = —2J 
and generates the TV-partite entangled W-state. (Let us 
mention in passing that the concurrence, which measures 
entanglement between any two qubits in the chain, has 
the value 2/N. In this sense the W state also carries bi- 
partite long-range entanglement.) This method does not 
require any dynamical control over the couplings in the 
spin chain. Only one external global parameter, namely 
the magnetic field B, has to be modified. Thus, our result 
opens new possibilities for the creation of multipartite 
entanglement in condensed matter physics. Moreover, if 
each qubit in the chain can be controlled separately then 
this transition can be achieved by local operations, i.e. 
by reducing the field locally at each on the N sites of 
the chain. This does not imply that entanglement can 
be created by local operations as qubits do interact with 
each other, but it opens new possibilities of manipulating 
multipartite quantum entanglement. 

We wish to thank M. Lewenstein and M. Palma for 
discussions. This work was supported in part by a 
grant from the Deutsche Forschungsgemeinschaft, the 
Cambridge-MIT Institute, A*Star Grant No. 012-104- 
0040 and the EU under project RESQ and QUPRODIS. 
We acknowledge support from the Benasque Center of 
Science during the Workshop on Quantum Information 
and Communication, 2003. 



[1] J. J. Bollinger et al, Phys. Rev. A 54, 4649 (1996). 

[2] A.K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

[3] C.H. Bennett et al, Phys. Rev. Lett. 70, 1895 (1993). 

[4] T. Yamamoto et al, Phys. Rev. A 66, 064301 (2002). 

[5] M. Eibl et al, Phys. Rev. Lett. 90, 200403 (2003). 

[6] S. Bose, Phys. Rev. Lett. 91, 207901 (2003), M. Chris- 
tandl et al, Phys. Rev. Lett. 92, 187902 (2004); C. Al- 
banese et al quant-ph/0405029, V. Subrahmanyam, Phys. 
Rev. A 69, 034304 (2004). 



[7] S.C. Benjamin and S. Bose, Phys. Rev. Lett. 90, 247901 

(2003) . 

[8] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 

416, 608 (2002). 
[9] F. Verstraete, M. Popp, and J.I. Cirac, Phys. Rev. Lett. 

92, 027901 (2004). 
[10] T. Osborne and M. Nielsen, Phys. Rev. A 66, 032110 

(2002) ; Y. Chen, P. Zanardi, Z.D. Wang, and F.C. Zhang, 
quant-ph/0407228; G. Vidal, J. I. Latorre, E. Rico, and 
A. Kitaev, Phys. Rev. Lett. 90, 22790 (2003); J. I. La- 
torre, E. Rico, and G. Vidal, Quant. Inf. and Comp. 4, 
48 (2004). 

[11] D. Gunlycke, S. Bose, V. M. Kendon, and V. Vedral, 
Phys. Rev. A 64, 042302 (2001); M. Asoudeh and V. 
Karimipour, quant-ph/0405060; X. Wang and P. Zanardi, 
Phys. Lett. A 301, 1 (2002); L. Zhou, H.S. Song, and 
Y.Q. Guo, Phys. Rev. A 68, 024301 (2003). 

[12] U. Dorner, P. Fedichev, D. Jaksch, M. Lewenstein, and 
P. Zoller, Phys. Rev. Lett. 91, 073601 (2003). 

[13] B.-Q. Jin and V.E. Korepin, J. Stat. Phys. 116, 79 

(2004) . 

[14] T. Roscilde, P. Verrucchi, A. Fubini, S. Haas, and V. 

Tognetti, Phys. Rev. Lett. 93, 167203 (2004). 
[15] X. Peng, J. Du, and D. Suter, Phys. Rev. A 71, 012307 

(2005) . 

[16] X. Wang, Phys. Rev. A 66, 044305 (2002). 
[17] P. Stelmachovic and V. Buzek, Phys. Rev. A 70, 032313 
(2004). 

[18] L.F. Santos, Phys. Rev. A 67, 062306 (2003). 

[19] D. Brufi and C. Macchiavello, Found. Phys. 33, 1617 

(2003) . 

[20] A. Acm, D. Brufi, M. Lewenstein and A. Sanpera, Phys. 
Rev. Lett. 87, 040401 (2001). 

[21] E. Lieb et al., Ann. Phys. (NY) 16, 407 (1961); S. Kat- 
sura, Phys. Rev. 127, 1508 (1962). 

[22] S. Sachdev, Quantum phase transitions, Cambridge Uni- 
versity Press (1999). 

[23] T. Kato, J. Phys. Soc. Jap. 5, 435 (1951). 

[24] V.N. Gorbachev, A.I. Trubilko, and A.A. Rodichkina, 
Phys. Lett. A 314, 267, (2003); J. Joo, Y.-J. Park, S. Oh, 
and J. Kim, New Jour, of Phys. 5, 136, (2003); Nguyen 
Ba An, Phys. Rev. A 69, 022315 (2004); D. Brufi et al., 
Phys. Rev. Lett. 93, 210501 (2004). 



